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Abstract 

Given a set 5 of 71 points in R^, and an integer k such that ^ fc < n, we show that a geometric 
graph with vertex set 5", at most n — 1 + k edges, maximum degree five, and dilation 0{n/{k + 1)) can be 
computed in time 0(n log n). For any k, we also construct planar n-point sets for which any geometric 
graph with n — 1 + k edges has dilation Q{n/{k + 1)); a slightly weaker statement holds if the points of 
S are required to be in convex position. 

1 Preliminaries and introduction 

A geometric network is an undirected graph whose vertices are points in M^. Geometric networks, especially 
geometric networks of points in the plane, arise in many applications. Road networks, railway networks, 
computer networks — any collection of objects that have some connections between them can be modeled 
as a geometric network. A natural and widely studied type of geometric network is the Euclidean network, 
where the weight of an edge is simply the Euclidean distance between its two endpoints. Such networks for 
points in M.^ form the topic of study of our paper. 

When designing a network for a given set S of points, several criteria have to be taken into account. In 
particular, in many applications it is important to ensure a short connection between every two points in 
S. For this it would be ideal to have a direct connection between every two points; the network would then 
be a complete graph. In most applications, however, this is unacceptable due to the high costs. Thus the 
question arises: is it possible to construct a network that guarantees a reasonably short connection between 
every two points while not using too many edges? This leads to the well-studied concept of spanners, which 
we define next. 

The weight of an edge e = {u, v) in a Euclidean network G ~ (S, E) on a set 5* of n points is the Euclidean 
distance between u and v, which we denote by d{u,v). The graph distance dG{u,v) between two vertices 
u,v £ S \s the length of a shortest path in G connecting uto v. The dilation (or stretch factor) of G, denoted 
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A(G'), is tlic maximum factor by wliicli the graph distance da differs from tlie Euclidean distance d, namely 



A(G) max — — . 

u,v<^s d[u,v) 

The network G is a t-spanner for S if A(G) ^ t. 

Spanners find applications in robotics, network topology design, distributed systems, design of parallel 
machines, and many other areas and have been a subject of considerable research [3j. Recently spanners 
found interesting practical applications in metric space searching [211 122j and broadcasting in communication 
networks [U [121 [H] ■ The problem of constructing spanners has received considerable attention from both 
an experimental perspective [TTl [M] and theoretical perspective — see the surveys [51 [131 Hi] and the book by 
Narasimhan and Smid [20j . 

The complete graph has dilation 1, which is optimal, but we already noted that the complete graph is 
generally too costly. The main challenge is therefore to design sparse networks that have small dilation. 
There are several possible measures of sparseness, for example the total weight of the edges or the maximum 
degree of a vertex. The measure that we will focus on is the number of edges. Thus the main question we 
study is this: Given a set 5 of n points in R^, what is the best dilation one can achieve with a network on 
S that has few edges? Notice that the edges of the network are allowed to cross or overlap. 

This question has already received ample attention. For example, there are several algorithms [21 \17\ 
[23l [26] that compute a (1 + £)-spanner for S, for any given constant e > 0. The number of edges in these 
spanners is 0{n). Although the number of edges is linear in n, it can still be rather large due to the hidden 
constants in the O-notation that depend on e and the dimension D. Das and Heffernan [5] showed how to 
compute in O(nlogn) time, for any constant e' > 0, a t-spanner with (1 + e')n edges and degree at most 
three, where t only depends on e' and D. 

Any spanner must have at least n — 1 edges, for otherwise the graph would not be connected, and the 
dilation would be infinite. In this paper, we are interested in the case when the number of edges is close 
to n — 1, not just linear in n. This leads us to define the quantity A(S', k): 

A(S,k):= min A(G). 

V(G) = S 
\E{G)\=n-l + k 

Thus A{S, k) is the minimum dilation one can achieve with a network on 5* that has n — 1 + k edges. 

Klein and Kutz [15] showed recently that given a set of points S, a dilation t and a number A: > 0, it is 
NP-hard to decide whether A(S', k) ^ t. Cheong et al. [1] showed that even the special case fc = is still 
NP-hard: given a set of points S and a real value t > 1, it is NP-hard to decide whether a spanning tree of 
S with dilation at most t exists. But for the special case when the tree is required to be a star, Eppstein 
and Wortman [10] gave an algorithm to compute minimum dilation stars of a set of points in the plane in 
O(nlogn) expected time. 

In this paper we study the worst-case behavior of the function A(S', k): what is the best dilation one can 
guarantee for any set S" of rt points if one is allowed to use n — \ + k edges? In other words, we study the 
quantity 

5{n, k) := sup A(S', fc). 

SCK° 
\S\=n 

For the special case when the set S is in and is required to be in convex positior0 we define the quantity 
5c{n,k) analogously. 

The result of Das and Heffernan [5] mentioned above implies that, for any constant e' > 0, S{n,e'n) is 
bounded by a constant. We are interested in what can be achieved for smaller values of fc, in particular for 
^ fc < n. 

In the above definitions wc have placed the perhaps unnecessary restriction that the graph may use no 
other vertices besides the points of S. We will also consider networks whose vertex sets are supersets of S. 



set of points is in convex position if they all lie on the boundary of their convex hull. 
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In particular, we define a Steiner tree on 5 as a tree T with 5* C V{T). The vertices in V{T) \ S are called 
Steiner points. Note that T may have any number of vertices, the only restriction is on its topology. 

Our results. We first show that any Steiner tree on a set S oi n equally spaced points on a circle has 
dilation at least n/Ti. We prove in a similar way that (5(n, 0) ^ f'^- ~ 1- We remark that Eppstein [9] gave a 
simpler proof of a similar bound for 5{'n, 0); we improve this bound by a constant factor. 

We then continue with the case < A: < n. Here we give an example of a set S of n points in the 
plane for which any network with n — 1 + fc edges has dilation at least —\n/{k + 1)J — 1, proving that 
5{n, k) ^ ^\n/{k + 1)J — 1. We also prove that for points in convex position the dilation can be almost as 
large as in the general case, namely Sci'^i ~ fl{n/ {{k + 1) logn)). 

Next we study upper bounds. We describe an O(nlogn) time algorithm that computes for a given set 
S oi n points in MP and a parameter ^ fc < n a network of at most n — 1 + k edges, maximum degree 
five, and dilation 0{n/{k + 1)). Combined with our lower bounds, this implies that S{n,k) = Q(n/{k + 1)). 
In particular, our bounds apply to the case k = o(n), which was left open by Das and Heffernan [5]. 
Notice that if A: ^ n, then we have 1 ^ S{n,k) ^ S{n,n — 1), and thus 6{n,k) = 8(1). This means that 
(5(71, k) = 6(1 + n/{k + 1)) holds for any /c ^ 0. 

Our lower bounds use rather special point sets and it may be the case that more 'regular' point sets 
admit networks of smaller dilation. Therefore we also study the special case when S* is a point set with 
so-called bounded spread. The spread of a set of n points S is the ratio between the longest and shortest 
pairwise distance in S. We show that any set S with spread s admits a network of at most n — 1 + k edges 
with dilation 0{s/{k + 1)-'^/^). This bound is asymptotically tight for s = 0{n/{k + 1)^^^/"^) (otherwise 
0{n/ [k + 1)) is a better bound). This also leads to tight bounds for the case when S* is a n^^^ x • • • x n^/^ 
grid. 

Notation and terminology. Hereafter S will always denote a set of points in MP . Whenever it causes 
no confusion we do not distinguish an edge e = (m, v) in the network under consideration and the line 
segment uv. 

2 Lower bounds 

In this section wc prove lower bounds on the dilation that can be achieved with n — \ -\- k edges, for 
^ fc < n, by constructing high-dilation point sets in K^. Of course, these lower bounds apply as well in 
higher dimensions. 

2.1 Steiner trees 

We first show a lower bound on the dilation of any Steiner tree for S. The lower bound for this case uses 
the set S" of n points pi,p2, . . . ,Pn spaced equally on the unit circle, as shown in Fig. [IJa), and is based on 
similar arguments as in Ebbers-Baumann et al. [J (Lemma 3). 

Theorem 1 For any n > 1, there is a set S of n points in convex position such that any Steiner tree on S 
has dilation at least . ; , > — . 

sin(7r/n) ir 

Proof. Consider the set S described above and illustrated in Fig. [Ija). Let o be the center of the circle, and 
let T be a Steiner tree on S. First, we assume that o does not lie on an edge of the tree. 

Let X and y be any two points and let 7 and 7' be two paths from x to y avoiding 0. We call 7 and 
7' (homotopy) equivalent if 7 can be deformed continuously into 7' without moving its endpoints or ever 
passing through the point o, that is, if 7 and 7' belong to the same homotopy class in the punctured plane 
M? \ {o}. 

Let 7i be the unique path in T from pi to Pi+i (where Pn+i ■= Pi)- We argue that there must be at least 
one index i for which 7; is not equivalent to the straight segment pipi+i, as illustrated in Fig.[T]^a). 
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(a) (b) 
Figure 1: (a) The homotopy class of the path 7,;. (b) lUustratmg the proof of Lemma [TJ 

We argue by contradiction. Let F be the closed loop formed as the concatenation of 71, . . . ,7„, and let 
F' be the closed loop formed as the concatenation of the straight segments piPi^i, for i = \ . . .n. If 7^ is 
equivalent to piPi^i, for all i, then F and F' are equivalent. We now observe that, since F' is a simple closed 
loop surrounding o, it cannot be contracted to a point in the punctured plane (formally, it has winding 
number 1 around o). On the other hand, F is contained in the tree T ^ o (viewed as a formal unio of its 
edges) and hence must be contractible in \ {o}. Hence F and F' cannot be equivalent, a contradiction. 

Consider now a path 7^ not equivalent to the segment PiPi+i- Then 7^ must "go around" o, and so 
its length is at least 2. The distance between pi and Pi+i, on the other hand, is 2sin(7r/n), implying the 
theorem. 

Now consider the case where o lies on an edge of T. Assume for a contradiction that there is a spanning 
tree T that has dilation 1/ sin(7r/n) — e for some e > 0. Let o' be a point not on T at distance e/100 from o. 
Then we can use the argument above to show that there are two consecutive points whose path in T must 
go around o' . By the choice of o' such a path must have dilation larger than 1/ sin(7r/n) — e, a contradiction. 
□ 



2.2 The case fc = 

If we require the tree to be a spanning tree without Steiner points, then the path 7; in the above proof must 
not only "go around" o, but must do so using points pj on the circle only. We can use this to improve the 
constant in Theorem [U as follows. Let Pi and Pi+i be a pair of consecutive points such that the path 7^ 
is not equivalent to the segment PiPi^i. Consider the loop formed by 7,^ and Pi+iPi- It consists of straight 
segments visiting some of the points of S. Let C be the convex hull of this loop. The point o does not lie 
outside C (otherwise, the loop would be contractible in the punctured plane) and so there exist three vertices 
vi, V2, and V3 of C such that o G AV1V2V3. Since the loop visits each of these three vertices once, its length 
is at least the perimeter of AV1V2V3, which is at least 4 by Lemma [T] below. Therefore, we have proven 

Corollary 1 For any Ji > 1, 

, , , 4-2sin(7r/n) 2?! 

2sm(7r/rij tt 

Lemma 1 Any triangle inscribed in a unit circle and containing the circle center has perimeter at least 4. 

Proof. Let o be the circle center, and let a, 6, and c be three points at distance one from o such that o is 
contained in the triangle Aabc. We will prove that the perimeter p{Aabc) is at least 4. 



^Notice that T may not be properly embedded in the piano, that is, the edges of T may cross or overlap. However, viewed 
not as a subset of the plane, but rather as an abstract simplicial complex, T is certainly simply connected and F is a closed 
curve contained in it and thus contractible, within T, to a point. Therefore it is also contractible in the punctured plane, as 
claimed. 
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We need the following definition: For a compact convex set C in the plane and ^ 6* < tt, let w(C, 9) 
denote the width of C in direction 0. More precisely, if is the line through the origin with normal 
vector (cos 6*, sin 6*), then w{C,9) is the length of the orthogonal projection of C to ig. The Cauchy-Crofton 
formula allows us to express the perimeter p{C) of a compact convex set C in the plane as p{C) = 

w{c, e)de. 

We apply this formula to Aabc, and consider its projection on the line £g (for a given 9). Wc choose a 
coordinate system where £g is horizontal. By mirroring and renaming the points a, b, c, we can assume that 
b is the leftmost point, a is the rightmost point, that o lies below the edge ab, and that c does not lie to 
the right of a. This immediately implies that w{oc,9) ^ w{ob,9). Consider now the point a' obtained by 
mirroring a at o. Since c must lie below the line aa' and not left of 6, we have w{oc, 9) w{oa' , 9) ~ w{oa, 9). 
This implies 3w{oc, 9) ^ w{oa, 9) + 'w{ob, 9) + w{oc, 9), and therefore 

w{Aabc,9) = w{oa,9) +w{ob,9) 

= w{oa, 9) + w(ob, 9) + w{oc, 9) — w{oc, 9) 
^ l{w{oa,9) +w{ob,9) +w{oc,9)). 

Integrating 9 from to tt and applying the Cauchy-Crofton formula gives p{Aabc) ^ |(p(oa) +p{ob) +p{oc)). 
Since oa, ob, and oc are segments of length 1, each has perimeter 2, and so we have p{Aabc) ^ | • 6 = 4. □ 



2.3 The general case 

We now turn to the general case, and wc consider graphs with n — I + k edges for < fc < n. 
Theorem 2 For any n and any k with < k < n, 



TT L/c + IJ 



Proof. Our example 5* consists oi k + I copies of the set used in Theorem [TJ More precisely, we choose sets 
Si, for 1 ^ i ^ k + 1, each consisting of at least [n/ {k + 1)J points. We place the points in Si equally spaced 
on a unit-radius circle with center at (2m, 0), as in Fig.O The set S is the union of S*!, . . . , Sk+i] we choose 
the cardinalities of the Si such that S contains n points. 







Figure 2: Illustrating the point set 5* constructed in the proof of Theorem [21 

Let G be a graph with vertex set S and n — 1 + fc edges. We call an edge of G short if its endpoints lie 
in the same set Si, and long otherwise. Since G is connected, there are at least k long edges, and therefore 
at most n — 1 short edges. Since X] l"^*! = '^j ^^^^ implies that there is a set Si such that the number of 
short edges with both endpoints in Si is at most \Si \ — 1. Let G' be the induced subgraph of Si. If G' is 
not connected, then Si contains two points connected in G using a long edge, and therefore with dilation at 
least n. If G" is connected, then it is a tree, and so the argument of Corollary [T] implies that its dilation is 
at least 



2-sin(7r/[?i/(A: + l)j; 



sin(7r/[n/(fc + l)J) 
This implies the claimed lower bound on the dilation of G. 



1 



- 1. 



□ 
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2.4 Points in convex position 



The point set of Theorem [T] is in convex position, but works as a lower bound only for fc = 0. In fact, by 
adding a single edge (the case k = 1) one can reduce the dilation to a constant. Now consider n points that 
lie on the boundary of a planar convex figure with aspect ratio at most p, that is, with the ratio of diameter 
to width at most p. It is not difficult to see that connecting the points along the boundary — hence, using 
n edges — leads to a graph with dilation Q{p). However, the following theorem shows that for large aspect 
ratio, one cannot do much better than in the general case. 

Theorem 3 For any n and any k with ^ k < n, 

Scin,k)^n(- V 

Proof. We present the proof for k — 0, but the construction can be generalized to hold for any A: > by 
using the same idea as in the proof of Theorem [51 placing fc + 1 copies of the construction along a horizontal 
line, with sufficient space between consecutive copies. 

For the case when fc = 0, set m = [n/4 — 1/2J and let o := (0,0). Consider the function f{i) = 
(1 + ^^^y ~ 1- Our construction consists of the \m + 2 ^ n points S with coordinates (/(«), 1), (/(i)j 
(— /(i), 1), and (— /(i), — 1), for i € {0, ... , m}, as shown in Fig. [3l The points of S lie on the boundary of 
a rectangle, so S is in convex position. (A slight perturbation of S would even give a set in strictly convex 
position, that is, a set where every point is extreme.) Consider a minimum-dilation tree T of S. As in the 
proof of Lemma [TJ we can now argue that there must exist two consecutive points p and q in S such that 
the path in T connecting them is not (homotopy) equivalent to the straight-line segment between them in 
the punctured plane K^\{o}. By symmetry, we can assume that p and q lie in the half-plane x ^ 0, and that 
p lies above the x-axis. 

The first case is p = (/(«), 1) and q = (/(i + 1), 1), where ^ i < m. The Euclidean distance between 
them is f{i + 1) - f{i) = ^-^{1 + f{i})), and the length of the shortest path in T is at least ^/2{l + f{i)). 
The two bounds imply that the dilation of T is at least v^j^ in this case. 

The second case is p = {f{m), 1) and q = {f{m), —1). Then the Euclidean distance between them is 2, 
and the length of the shortest path in T between them is at least 2/(to), thus T has dilation /(m) in this 
case. It remains to bound f{m), which can be done by using the inequality (1 + t/m)"^ ^ e*(l — t^ /m), 
which holds for |t| ^ m [HI Proposition B.3]. We obtain 

„ , . / In TO \ / In^ TO \ ^ , X 

/(m)= IH -l^m 1 -l = f7(m), 

V TO / V TO / 

which concludes the proof of the theorem. H 

(-./(Oa) •'^ = (/(*), 1) 



y = 1 • 




y = — 1 • • • ••••••••••• • • • 

(-/(*),-!) ■ (/(*),-!) 

Figure 3: Illustrating the lower bound construction used in the proof of Theorem [31 with n = 34. 
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3 A constructive upper bound 



In this section wc show an upper bound on the dilation achievable with k extra edges. We make use of the 
fact (observed by Eppstein [5]) that a minimum spanning tree has dilation at most n — 1. We will use the 
following lemma. 

Lemma 2 Let S be a set of n points in M.^ , and let T be a minimum spanning tree of S. Then T has 
dilation at most n — 1. Hence A{S,0) ^ A(T) ^ n — 1 and, as this holds for all S with \S\ = n, we have 
(5(n,0) 

Proof. Let p,q Cz S and consider the path 7 connecting p and q in T. Since T is a minimum spanning tree, 
any edge in 7 has length at most d{j},q). Since 7 consists of at most n — 1 edges, the dilation of 7 is at 
most n — 1. □ 

3.1 The planar case, D = 2 

The following algorithm builds a spanner with at most n — 1 + fc edges: 



Algorithm 1 SparseSpanner(S', fc) 

Input A set S" of n points in the plane and an integer fc € {0, . . . , 2n — 5}. 

Output A graph G = {S, E) with dilation 0{n/ (fc + 1)) and at most n — 1 + fc edges. 

1: Compute a Delaunay triangulation of S. 

2: Compute a minimum spanning tree T of S. 

3: if fc = then 

4: return T. 

5: Let [(fc + 5)/2j. 

6; Compute m disjoint subtrees of T, each containing 0(n/m) points, by removing m — 1 edges. 
7: E^9. 

8; for each subtree T' do 
9: add the edges of T' to E. 
10: for each pair of subtrees T' and T" do 

11: if there is a Delaunay edge {p, q) with p £ T', q £ T" then 
12: add the shortest such edge {p, q) to E. 
13: return G = {S,E). 



We first prove the correctness of the algorithm. 

Lemma 3 Algorithm SparseSpanner returns a graph G with at most n — l + k edges and dilation bounded 
byO{n/{k + l)). 

Proof. Lemma [2] shows that our algorithm is correct if fc = 0, so from now on we assume that fc 1, and 
thus m ^ 3. The output graph G is a subset of a Delaunay triangulation of S and is therefore a planar 
graph, see Fig. |4^. Consider now the graph G' obtained from G by contracting each subtree T' created in 
step 6 to a single node. Since G is planar and since two subtrees are connected with at most one edge, G' is 
a planar graph with m ^ 3 vertices, without loops or multiple edges, and so it has at most 3to — 6 edges. 
The total number of edges in the output graph is therefore at most 

(71 - 1 - (m - 1)) + (3m -6)=n + 2m-6^n + fc-l. 

Next we prove the dilation bound. Consider two points x,y G S. Let x = xq,xi, . . . ,Xj = y he a. shortest 
path from a; to y in the Delaunay graph. The dilation of this path is bounded by 27r/(3 cos(7r/6)) = 0(1) [l4] . 
We claim that each edge XiXi+i can be replaced by a path in G with dilation 0{n/k). The concatenation of 
these paths yields a path from x to y with dilation 0{n/k), proving the lemma. 
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If Xi and Xi^i fall into the same subtree T', then Lemma [5] implies a dilation of 0{n/m) ~ 0{n/k). 

It remains to consider the case Xi S T', Xi+i S T", where T' and T" are distinct subtrees of T, see Fig.|4}3. 
Let (p, (7) be the edge with p G T' , q G T" inserted in step 12 — such an edge exists, since there is at least one 
Delaunay edge between T' and T", namely {xi,Xi+i). Let fi' be a path from Xi to p in T', and let fi" be a 
path from g to x^+i in T". Both paths have 0{n/k) edges. By construction, we have d{j),q) ^ d{xi^Xi^i). 
We claim that every edge e on /x' and /i" has length at most d{xi^ Xi+i). Indeed, if e has length larger than 
Xi+i), then the original tree T (of which T' and T" arc subtrees) cannot be a minimum spanning tree 
of S: we can remove e from T and insert either [xi, Xi+i) or (p, to obtain a better spanning tree. Therefore 
the concatenation of /x', the edge (p, g), and /x" is a path of 0(n/k) edges, each of length at most d{xi, Xi+i), 
and so the dilation of this path is 0{n/k). H 



Figure 4: (a) The solid edges connect vertices within a subtree, while the dotted edges connect different 
subtrees, (b) Illustrating the proof of Lemma [31 

Theorem 4 Given a set S of n points in the plane and an integer ^ fc ^ 2n — 5, a graph G with vertex 
set S*, n — 1 + fc edges, and dilation 0(n/{k + 1)) can he constructed in O(nlogn) time. 

Proof. We use algorithm SparseSpanner. Its correctness has been proven in Lemma [31 Steps 1 and 2 
can be implemented in O(nlogn) time [9]. Step 6 can be implemented in linear time as follows. Orient T 
by choosing any node of degree at most five as root node. Since the degree of a minimum spanning tree is 
at most six I27|. every node in T has at most five children. We will partition T into disjoint subtrees by 
removing edges, such that each subtree has size at least n/m and at most 5(n/m) — 4. The only exception 
is the top subtree — that is, the subtree containing the root r of T — which has size at most 5{n/m) — 4 but 
could be smaller than n/m. To obtain this partition, we first partition the trees rooted at the at most five 
children of r recursively. All subtrees wc get in this manner will have sizes between n/m and 5{n/m) — 4, 
except possibly for the subtrees containing the children of r which could be too small. The subtrees that are 
too small are grouped together with r into one new subtree of size at most b{n/m) — 4. Since all but one 
subtree have size at least n/m, the number of subtrees is at most m. To obtain exactly m subtrees, we can 
simply remove some more edges arbitrarily. To implement steps 10-12, we scan the edges of the Delaunay 
triangulation, and keep the shortest edge connecting each pair of subtrees. This can be done in 0{nlogn) 
time. H 

3.2 The higher-dimensional case, D ^ 2 

The algorithm of Section 13.11 uses the fact that the Delaunay triangulation has constant dilation, a linear 
number of edges, and can be computed efficiently. These properties do not generalize to dimensions D ^ 3, 
and so we have to turn to another bounded-degree spanner of our point set S. Since the minimum spanning 
tree is unlikely to be computable in near- linear time in dimension D ^ 3, wc will use a minimum spanning 
tree of this spanner instead of a minimum spanning tree of the point set itself. These two ideas will allow 
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us to generalize our result to any constant dimension D ^ 2. Moreover, we show that this can be achieved 
by a graph with degree at most five. 

Properties of the minimum spanning tree of a spanner. Let i ^ 1 be a real number, let G be an 
arbitrary t-spanner for S, and let T be a minimum spanning tree of G. In the following three lemmas, we 
prove that T has "approximately" the same properties as an exact minimum spanning tree of the point set S. 

Lemma 4 Let p and q be two points of S . Then every edge on the path in T between p and q has length at 
most t ■ d{p, q). 

Proof. Let (x, y) be an arbitrary edge on the path in T from p to q. For a contradiction, assume that 
d{x,y) > t ■ d{p,q). 

Let T' be the graph obtained by removing the edge (x, y) from T. It consists of two components, one 
containing p, the other containing q. Now, since G is a i-spanner, it contains a path 7 between p and q of 
length at most t ■ d{p, q). The union of T' with the path 7 is a spanning graph of G, and its weight is less 
than the weight of T, a contradiction. H 

Lemma 5 T is a t{n — \)-spanner for S . 

Proof. Let p and q be distinct points of 5, and let 7 be the path in T between p and q. By Lemma |H each 
edge of 7 has length at most t ■ d{p, q). Since 7 contains at most n — 1 edges, it follows that the length of 7 
is at most t{n — 1) • d{p, q). H 

Lemma 6 Let m be an integer with 1 ^ to ^ n — 1, and let T' and T" be two vertex-disjoint subtrees of T , 
each consisting of at most m vertices. Let p be a vertex of T' , let q be a vertex of T" , and let 7 be the path 
in T between p and q. If x is a vertex of T' that is on the subpath of ^ within T' , and y is a vertex of T" 
that is on the subpath of j within T" , then 

d{x,y) ^ (2t(m-l) + l)-d{p,q). 

Proof. Let 7' be the subpath of 7 between p and x. By Lemma HI each edge of 7' C 7 has length at most 
t ■ d{j), q). Since 7' contains at most to — 1 edges, it follows that this path has length at most t(rn — 1) ■ c?(p, q). 
On the other hand, since 7' is a path between p and x, its length is at least d{p,x). Thus, we have 
d{p, x) ^ t{m — 1) • d(j), q). A symmetric argument gives d{q, y) ^ t(m — 1) • d{p, q). Therefore, we have 

d{x, y) ^ d{x,p) + d{p, q) + d{q, y) 

^ t{m - 1) • d{p, q) + d{p, q) + t{m - 1) • d{p, q), 

completing the proof of the lemma. H 

A graph with at most n — 1 + k edges and dilation 0{n/{k + 1)). Fix a constant t > 1. Let c be a 
constant such that for any set X of 2m points, the t-spanner of Das and Heffernan [S] has at most cm edges. 
Let G be a i-spanncr for 5* whose degree is bounded by a constant. Since then the minimum spanning tree 
T of G has bounded degree as well, we can use the algorithm described in the proof of Theorem [1] (step 6) 
to remove to ~ lk/{c — 1)J edges from T and obtain vertex-disjoint subtrees Tq, Ti, . . . , T^, each containing 
0{n/{k + 1)) vertices. The vertex sets of these subtrees form a partition of S. Let X be the set of endpoints 
of the m edges that are removed from T. The size of X is at most 2m. 
We define G' to be the graph with vertex set S that is the union of 

1. the trees Tq, Ti, . . . , Tm, and 

2. Das and Heffernan's i-spanner G" for the set X (G" is empty if m = 0). 

We first observe that the number of edges of G' is bounded from above by n — 1 — to + cm ^ n — 1 + k. 
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Lemma 7 The graph G' has dilation 0{n/{k + 1)). 

Proof. The statement follows from Lemma [S] if m ~ 0. Let m > 0, and let p and q be two distinct points 
of S. Let i and j be the indices such that p is a vertex of the subtree Ti and g is a vertex of the subtree Tj. 

First assume that i = j. Let 7 be the path in Ti between p and q. Then, 7 is a path in G' . By LemmaHl 
each edge on 7 has length at most t ■ d{p, q). Since Ti contains 0{n/k) vertices, the number of edges on 7 is 
0{n/k). Therefore, since t is a constant, the length of 7 is 0{n/k) ■ d{p, q). 

Now assume that i ^ j. Let 7 be the path in T between p and q. Let (x, x') be the edge of 7 for which x 
is a vertex of Ti, but x' is not a vertex of Ti. Similarly, let (y, y') be the edge of 7 for which y is a vertex of 
Tj, but y' is not a vertex of Tj. Then, both {x, x') and [y, y') are edges of T that have been removed when 
the subtrees were constructed. Hence, x and y are both contained in X and, therefore, are vertices of G" . 
Let 7j be the path in Tj between p and x, let 7^,^ be a shortest path in G" between x and y, and let 7^ be 
the path in Tj between y and g. The concatenation 7' of 7^, ^xy, and 7^ is a path in G" between p and q. 

Since both 7^ and 7^ are subpaths of 7, it follows from Lemma |4] that each edge on 7^ and 7^ has length 
at most t ■ d{p, q). Since Ti and Tj contain 0{n/k) vertices, it follows that the sum of the lengths of 7^ and 
7j is 0{n/k) ■ d{p,q). The length of jxy is at most t ■ d{x,y) which, by Lemma [6l is also 0(n/k) ■ d{p,q). 
Thus, the length of 7' is 0(n/fc) • d{p, q). □ 

Now let G be the i-spanner of Das and Heffernan [S]. This spanner can be computed in 0(?ilogn) 
time, and each vertex has degree at most three. Given G, its minimum spanning tree T and the subtrees 
Tq,Ti, . . . ,Tm can be computed in O(nlogn) time. Finally, G" can be computed in O(TOlogm) = 0{n\ogn) 
time, and each vertex has degree at most three. Thus G" has dilation 0{n/{k + 1)), it contains at most 
n — 1 + k edges, and it can be computed in 0{n log n) time. 

We analyze the degree of G' : Consider any vertex p oi G' . li p ^ X, then the degree of p in G' is equal 
to the degree of p in T, which is at most three. Assume that p E X. The graph G" contains at most three 
edges that are incident to p. Similarly, the tree T contains at most three edges that are incident to p, but, 
since p G X, at least one of these three edges is not an edge of G". Therefore, the degree of p in G" is at 
most five. Thus, each vertex of G' has degree at most five. 

Thus, we have proved the following result: 

Theorem 5 Given a set S of n points in and an integer k with ^ k < n, a graph G with vertex set S , 
n — 1 + k edges, degree at most five, and dilation 0{n/{k + 1)) can be constructed in O(nlogn) time. 

4 Bounded Spread 

In this section we consider the case when the set of input points has bounded spread. The spread of a set 
of n points S, denoted s{S), is the ratio between the longest and shortest pairwise distance in S. In we 
haves(5) =17(ni/-°). 
We define the function 

S{n, s, k) := sup{A(S', fc) | 5 C M^, |5| = n, s{S) s}, 

that is, the best dilation one can guarantee for any set 5 of n points with spread s if one is allowed to use 
n — 1 + k edges. 

Theorem 6 For any n, any s, and any k with ^ k < n, 

S{n,s,k) = 0{s/{k + iy/^). 

Proof. Assume without loss of generality that the smallest interpoint distance in 5 is 1. 

The case fc 2 is nearly trivial: we pick an arbitrary point u Cz S and connect all other points to u. Since 
the shortest interpoint distance is 1, the dilation is at most 2s. 

For the case k > 2, let B be an axis-parallel cube with side length s containing S, and let m := [(fc— 1)^/-'^J . 
We partition B into small cubes with side length s/m. Let us call each small cube a cell. Let X C S he 
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a set of representative points, that is, one point of S taken from each non-empty cell. Let G' be a t-spanner 
for X with 2\X\ edges and constant t [5]. We obtain our final graph G from G' by connecting each point 
p G 5* to the representative point of the cell containing p. 

Since \X\ ^ ^ fc — 1, the total number of edges of G is at most n — \X\ + 2\X\ ^ n — 1 + fc. 

It remains to prove that the dilation of G is 0{s/m) for every two points p,q € S. If p and q lie in the 
same cell, then the bound follows immediately. If p and q lie in different cells, then let p' and q' be their 
representative points. The length of the shortest path between p and q in G is then 

dciP, q) = dcip.p') + dG{p\q') + doiq^q) 

< d{p,p ) + t ■ d{p,q ) + d{q,q) 

< d{p,p') + t ■ {d{p',p) + dip, q) + d{q, q')) + d{q' , q) 
= (1 + t){d{p,p') + d{q', q))+t- dip, q) 

^ {l + t)2VD—+t-d{p,q) 
m 

< {{l + t)2y/D— +t) ■d{p,q). 



The last inequality uses d{p, q) ^ I. H 
To prove a corresponding lower bound, we need a lemma in the spirit of Theorem [1] 

Lemma 8 Let Q be a two-dimensional square of side length 2a in MP , let S be a set of points on the 
boundary of Q, and let T be a Steiner tree for S (with vertices in M.^ ). Then there are two consecutive points 
p, q € S such that drip, q) ^ 2a. 

Proof. Let T' be the orthogonal projection of T onto the plane spanned by Q. Edges of T may cross 
or overlap in the projection, vertices may even coincide, but viewed as a formal union of its edges. T' is 
now a Steiner tree for S. The argument made in the proof of Theorem [1] therefore goes through; there are 
consecutive points p, q in S such that the path from p to q in T' goes around the center of Q, and therefore 
its length is at least 2a. H 

Theorem 7 For any integer r ^ 4, any integer m with 1 ^ m ^ r/4, and any integer n with 2rm^^^ ^ 
n ^ r^ , there is a subset S of n points of the r x r x ■ ■ ■ x r integer grid {0, 1, . . . , r — 1}^ such that 

2m 

Proof. Let a := [r/(4m)J . We define small subgrids as follows: For an integer vector x = (ii, ^2, ■ • ■ , io) G 
{0, ...,m — 1}^, let the cell G{x) := {0,1,..., 4ct — 1}-^ + Aax. For cell G{x), we define its square 
Q{x) := {{a, . . . ,3(t}^ x {2(t}^~^) + 4:ax, as illustrated in Fig. O Note that cells are D-dimensional and 
contain 4'^ct^ points each, squares are two-dimensional and contain {2a + 1)^ points each. We are in fact 
only interested in the boundary points of each square, there are Scr such boundary points in each square, 
and therefore 8m^a such points in total. Our set S contains all these square boundary points, and then an 
arbitrary subset of the remaining grid points {0, . . . , r — 1}^. 

It remains to prove that any network G with n — 1 + {m^ — 1) edges on 5* has dilation at least 2a. Let 
us call an edge e of G short if both endpoints of e are in the same cell, and long otherwise. Let T be any 
spanning tree of G. We color the — 1 edges of G that are not in T red. Since there are at most — 1 
short red edges in G and cells, there must be a cell C(x) without a short red edge. We apply Lemma [5] 
to the boundary points S' of the square Q{x) inside C{x). The tree T is a Steiner tree for S' , and so there 
are consecutive points p, q £ S' such that drip, q) ^ 2a. If daip, q) = drip, q), then this implies A(G) ^ 2a. 
Otherwise, the shortest path from p to g in G must use a red edge. Since C{x) contains no short red edge, 
the path must therefore pass through a point outside C{x), which implies that it has length at least 2a, and 
so A(G) ^ 2a. □ 

Combining Theorems [6] and [7l we obtain an asymptotically sharp bound for the dilation. 
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(0,4iT- 1,4(T- 1) +4ax 
(4ct-1,0,4<t-1)+4cti: 



2a + 4ai, 




(4(t - 1, 4<T - 1,4(7 - 1) +4(ti: 



^ (4ct - 1,4<t - 1,0) +4CTa: 



(0,0,0) +4(Ti: 



(4cr - 1,0,0) +4(ja: 



Figure 5: The cube illustrates the cell of an integer vector a; in 3 dimensions and the set of solid black points 
illustrates the square of the cell. 

Corollary 2 For any n, any k with ^ k < n, and any s with s ~ Q{n^/^) and s = 0{n/{k + 1)^^^/^), 
we have 

6in,s,k) = e(s/(fc + 1)1/-°). 

Proof. The upper bound is Theorem [6l For the lower bound, set r := [s/a/D] and m \{k + 1)^/^], 
apply Theorem [7] and observe that the resulting grid set has spread at most s. H 

The "most regular" point set one might imagine is the regular grid. It turns out that the spread-based 
lower bound remains true even for this set. 

Corollary 3 Let Q be a set of n = r^ points forming a r x r x ■ ■ ■ x r-grid, and let ^ k < n. Then 

A(5,fc) = e(r/(fc + l)i/^). 

Proof. The upper bound follows from the fact that Q has spread s = 0(r) (better constants can be obtained 
by a direct construction) . The lower bound follows from Theorem [71 by observing that if you set n to , 
then the construction results in S = Q. H 



5 Open problems 

We have shown that for any n-point set S in and any parameter ^ k < n, there is a graph G with 
vertex set 5, n — 1 + fc edges, degree at most five, and dilation 0{n/ {k + 1)). We also proved a lower bound of 
n{n/{k + 1)) on the maximum dilation of such a graph. An interesting open problem is whether the degree 
can be reduced to four, or even three. 

The constant in our lower bound for 5{n/k) is 2/-K, but we have only proven asymptotically matching 
upper bounds. Even for the case fc = 0, the upper bound is n — 1 while the lower bound is only 27i/7r — 1, 
and it would be interesting to establish the right constant. For fc > the discrepancy between upper and 
lower bounds is even larger. 

Minimum-dilation graphs are not well understood yet. As mentioned in the introduction, it is NP-hard 
to decide whether there is a t-spanner of a point set with at most n — 1 -I- fc edges, even in the case when 
the spanner is a tree [?]. However, if the spanner is restricted to be a star then the minimum dilation graph 
can be computed in polynomial time |10| . Are there other restrictions on the spanner or the point set that 
allow for efhcient algorithms? 

Given that the general problem is NP-hard, it would also be interesting to look for algorithms that 
approximate the best possible dilation (instead of giving only a guarantee in terms of n and fc, as we 
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do). We are not aware of any result showing how to approximate the minimum-dilation spanning tree 
with approximation factor o{n). For general unweighted graphs, a O(logn) approximation is possible [8]. 
Another result in this direction is by Knauer and Mulzer |16j . who described an algorithm that computes a 
triangulation whose dilation is within a factor of 1 + 0(1/ y/n) of the optimum. 
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